Juxtaposition of Spin Freezing and Long Range Order in a Series of Geometrically 
Frustrated Antiferromagnetic Gadolinium Garnets 
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Specific heat measurements in zero field are presented on a homologous series of geometrically 
frustrated, antiferromagnetic, Heisenberg garnet systems. Measurements of Gd3Ga50i2, grown 
with isotopically pure Gd, agree well with previous results on samples with naturally abundant 
Gd, showing no ordering features. In contrast, samples of Gd3Te2Li30i2 and Gd3Al50i2 are found 
to exhibit clear ordering transitions at 243 mK and 175 mK respectively. The effects of low level 
disorder are studied through dilution of Gd''^ with non-magnetic Y''^ in Gd3Te2Li30i2. We discuss 
possible explanations for such diverse behavior in very similar systems. 

PACS numbers: 75.50.Lk, 75.50.Ee, 75.40.Cx 



Magnetic systems of spins residing on the sites of a 
lattice of corner-sharing simplexes, such as triangles or 
tetrahedra, coupled via a nearest-neighbor (n.n.) an- 
tiferromagnetic Heisenberg exchange Hamiltonian, T-Lq, 
are highly frustrated. The spins are unable to minimize 
their energy pair by pair and, at the classical level, such a 
model results in a sort of "spin liquid" state without long 
range order (LRO) and with zero net magnetic moment, 
Ms, on each simplex Perturbations to Ho, 'H' , such 
as exchange beyond n.n. and dipolar interactions, per- 
haps assisted by thermal and/or quantum fluctuations, 
are expected to lift the classical ground state degeneracy 
and drive the system into a state of LRO. However, the 
ground state remains extremely fragile against quenched 
disorder and as a result, rather than developing LRO, the 
combination of high frustration and weak random disor- 
der can in principle 0, Q and in practice, for example 
in the pyrochlore Y2M02O7 (YMO) 0, [Hj, cause a sys- 
tem to exhibit a spin glass transition. One might quickly 
conclude that the glassy physics of the highly frustrated 
magnet (HFM) Gd3Ga50i2 (GGG), a hyperkagome net- 
work of corner-sharing triangles, is the result of a similar 
situation. However, the glassy behavior of YMO seems 
largely conventional 0] whereas GGG is known to dis- 
play highly exotic properties that incorporate unconven- 
tional spin glass behavior Q coexisting with extended 
short range order and low temperature spin fluctua- 
tions [i. 

Thus one is led to consider the much more exciting pos- 
sibility that the presence of several interaction terms of 
different spatial range and anisotropic nature that make 
up Ho and %' may, through their competition, lead to a 
complex energy landscape causing glassybehavior that 
is not induced by quenched randomness [9[ . Such a phe- 



nomenon parallels a theoretical description of structural 
glasses that builds on the notion of locally preferred 
structure in a liquid that is frustrated at large length 
scale due to its inability to suitably flll space, hence in- 
hibiting the formation of a crystalline state and leading 
to a glass transition [l^l . Could a HFM system, with its 
locally satisflcd Ms = simplexes, display a disorder-free 
freezing akin to the glass transition when subject to the 
competition of the terms that constitute T-L'l This may 
be well realized in a system like GGG where the domi- 
nant local Heisenberg exchange is in competition with the 
also important long range dipolar interaction [ll. 12|. 

The amount of disorder in GGG is low, consisting of 
up to a possible 1-2% off-stoichiometry of Gd''"*' magnetic 
ions on otherwise non-magnetic Ga'^"'' sites, inherent in 
crystal growth [13|. There has so far been no theoreti- 



cal explanation for the glassiness seen in GGG, thus it 
is crucial to ascertain whether the behavior in GGG is 
an intrinsic property of Gd^+ ions on the garnet lattice 
arising from competing antiferromagnetic short range ex- 
change and long range dipolar interactions or whether a 
generic 1-2% random disorder is responsible for the glass 
transition. In this Letter, we investigate this question 
using specific heat measurements to study a homologous 
scries of three Gd garnet materials: the aforementioned 
GGG as wefl as GdgTeaLigOia (GTLG) and GdgAlsOiz 
(GAG). Additionally, we have tested the effects of ran- 
dom chemical substitution of non-magnetic Y'^"'' for Gd'^"'' 
in GTLG to investigate in a controlled manner a level of 
disorder comparable to or higher than that of GGG. Our 
results show that spin glass behavior resulting from a 
generic sensitivity to disorder is too simple a picture for 
this particular series of materials. 

Wc have measured a single crystal of GGG, made from 
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FIG. 1: (color online) Specific heat (C) and C/T of isotopi- 
cally pure GGG measured with coarse temperature resolution. 
Results match well with previous work [^. The temperature 
region marked with a red box has been remeasured with a high 
temperature resolution shown in the inset and again does not 
find evidence of an ordering transition. 



isotopically pure ^^"Gd, grown from the same powder 
used for neutron scattering measurements 0, [13. Prepa- 
ration of this sample has been described in Ref. [l5| . The 
samples of GTLG and GAG were polycrystalline samples 
made from naturally abundant Gd. GTLG samples were 
made by solid-state reaction in air of Te and Gd (Y) ox- 
ides and Li carbonate, pressed into pellets, at 850 C for 
10 hours, then rcgrinding and firing a second time. The 
polycrystalline sample of GAG was made using the sol-gel 
method It was heated to 1350 C for one hour, then 
removed from the furnace and rapidly quenched to room 
temperature, to minimize the formation of perovskite- 
phase GdAlOa Powder x-ray diffraction spectra in 
both cases could be indexed to space group laSd. No 
impurity peaks were visible at the 1% level in GTLG; 
in GAG, there were two impurity peaks attributed to 
GdAlOa, at the 3% level. 

Specific heat, C, measurements were performed at 
dilution refrigerator temperatures, using the quasi- 
adiabatic method as in Refs. 0, [lali with thermometer 
and heater fixed directly to the samples. 6 /xm diameter, 
^ 1 cm long superconducting leads were used to provide 
excellent thermal isolation. A long time constant of re- 
laxation (r > 1 hour) ensured that internal temperature 
gradients due to poor thermal conductivity (of particu- 
lar concern in the powder samples) were not a significant 
source of systematic error. 



An initial measurement of C of GGG was made over 
a large temperature range (from 80 mK to 930 mK) us- 
ing temperature steps of roughly 5 mK below 200 mK 
and temperature steps of 10 mK above 200 mK. Results 
of our specific heat measurements on GGG, shown in 
Fig. [U agree remarkably well with the previous specific 
heat measurement (on a single crystal containing natu- 
rally abundant Gd) of Schiffer et al. Q . We find a broad 
feature with a maximum at around 800 mK. This fea- 
ture seems to drop out at lower temperature roughly as 
2^0.8 around 125 mK, at which point there is a max- 
imum in C/T and C{T) becomes steeper. The specific 
heat measurement of Dunsiger et al. on a naturally abun- 
dant Gd, powder sample is very similar, though the peak 
in C/T is found to be more pronounced Q. This result 
suggests that the isotopically pure sample measured here 
exhibits the same physics as do naturally abundant Gd 
containing samples and reconfirms the absence of a sharp 
ordering feature in C that would indicate a transition to 
LRO. 

Other experiments by Schiffer et al. on GGG have sug- 
gested the material to be a spin glass, albeit a somewhat 
unconventional one 0]. The spin glass interpretation is 
largely based on glassy relaxation in the ac susceptibil- 
ity and a sharp peak in the nonlinear susceptibility {x^) 
at around 180 mK Q. However, there is also a broader 
peak in at 450 mK and there is no corresponding max- 
imum in the specific heat (as seen in Fig. [T]) but just a 
maximum in C/T at around 120 mK. 

There is also a collection of experiments that pro- 
vide evidence against a spin glass transition. ^SR and 
Mossbauer experiments, for example, show significant 
persistent spin dynamics (PSDs) down to temperatures 
as low as 25 mK [8, 2^, 2l|. Most curiously, neutron scat- 
tering experiments on GGG, possible only with samples 
made from isotopically pure ^^"Gd that does not absorb 
neutrons, have shown sharp diffraction peaks developing 
below 140 mK Though these peaks are not sharp 

enough to imply true long range order (LRO), they sug- 
gest magnetic ordering with a correlation length of at 
least 100 A. This has been suggested to be a type of 
mixed spin liquid/solid state [7|. Bulk dynamical mea- 
surements have also been interpreted as a signature of an 
ordering transition [2^ . Recent inelastic neutron scatter- 
ing measurements have revealed three gapped dispersion- 
less excitations are found, two of which can be attributed 
to dimerized antiferromagnetic order [23I . 
In a theoretical study 



111 . [12], aiming to relate the 



sharp features in the GGG neutron scattering pattern Q 
with a long-range ordered state, the spin-spin correla- 
tions were treated using mean-field theory and using 
an Ewald summation method to handle the important 
dipolar interactions. Tuning the second (J2) and third 
(J3) n.n. couplings and simulating the powder neutron 



diffraction signal, Ref. [ll| was able to find excellent 



agreement with experiment, thereby determining optimal 
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FIG. 2: (color online) Specific heat of GTLG showing a 
sharp first-order phase transition. The specific heat drops 
out quicker than below Tc- Inset (a) shows the entropy 
as a function of T compared with the total _R In 8 entropy in 
the system. Inset (b) shows CT^ as a function of and a 
linear fit suggestive of spin waves with a gap of A = 0.62 K. 



interaction strengths. Most interestingly, it was found 
that the system exhibits a quasi-degeneracy critical (soft) 
modes, suggesting that this enhances fluctuations and 
makes it very sensitive to small amounts of disorder [l^l . 

This apparent development of rather long range or- 
der and its seeming consistency with theory [ll|, [l2j 
makes the lack of an ordering transition in specific heat 
measurements very paradoxical indeed. In order to ad- 
dress this problem, we have performed an additional mea- 
surement using a much higher temperature resolution of 
1 mK over the temperature range 130 mK to 230 mK to 
search for small or narrow features near where the neu- 
tron scattering peaks were discovered that might have 
been previously missed. This choice of resolution is based 
on scaling results of well characterized antiferromagnets 
to a transition temperature of 140 mK, suggesting one 
might expect a peak in C with a width of only several 
mK. However, the high-resolution scan of C, shown in 
the inset of Fig. [TJ also does not reveal any anomalies 
that might be interpreted as an ordering transition. 

In stark contrast, GTLG displays a sharp transition at 
243 mK as shown in Fig. [5] This is close to the tempera- 
ture 250 mK) where a feature was previously observed 



in the magnetic susceptibility of GTLG [2J| . The transi- 
tion is clearly first-order, exhibiting a much sharper peak 
in C than could be expected for a continuous phase tran- 
sition. With a Curie- Weiss temperature 9c w — —2.7 
in GTLG 2J1 , this ordering temperature gives a frustra- 
tion index / = 9cw/Tc = 11. Below the transition, C 
drops out faster than suggesting that it is exponen- 
tial and that there are gapped spin wave excitations, as 
would be expected from LRO with a strong dipolar in- 
teraction [3, HH . Gapped spin waves should result in a 
low temperature behavior C oc r~^e~^/^, thus we plot 
CT2 vs. 1/T in inset (b) of Fig. H The resulting lin- 
ear fit gives A ~ 0.62 K. Above the transition, there 
is a broad feature centered around roughly 1.0 K, sim- 
ilar to the broad feature in GGG centered at ~ 0.8 K 
(see Fig. [T]) and is likely a signature of the development 
of short range correlations. Obtaining the entropy (S) 
from a numerical integral oi C/T shows that only about 
14% of the total i? In 8 entropy in the system is accounted 
for by the transition, as shown in inset (a) of Fig. [51 

The third material studied here, GAG, may represent 
a point in between GGG and GTLG with a smaller and 
broader transition at a lower temperature of 175 mK 
(shown in Fig. ^. With 0cw ^ -3.0 K [l^, it is more 
antiferromagnetic than the other two garnets, but seem- 
ingly more frustrated than GTLG, with a frustration 
index / = 17. Otherwise, it shows similar features to 
GTLG, with a broad maximum centered around 1 K and 
a steeply dropping specific heat at lower temperatures. 
As in GTLG, the transition in GAG accounts for only a 
small percentage of the total R In 8 entropy in the system. 
The smaller and broader transition may be the result of 
the small level of GdAlOa impurities. 

It seems likely that GTLG does not have the same off- 
stoichiometry that is found in GGG, making it a "cleaner 
sample" ; Gd'^^ ions are unlikely to take the place of Li+ 
or Te^+ ions as easily as Ga^+ ions. Thus, we have added 
a low level of disorder to GTLG through dilution of Gd'^+ 
ions with non-magnetic Y'^^ ions. This "dirty" GTLG 
sample shows a significantly broadened peak but no sig- 
nificant change in the peak temperature, as shown in 
Fig. [3l At temperatures well below and well above the 
transition, C of the pure and diluted samples match, and 
the transition accounts for the same amount of entropy in 
both systems, suggesting that the ground state ordering 
is not appreciably altered by the 2% impurity. 

In conclusion, despite their similar models, we find 
GTLG and GAG to in fact behave entirely differently 
from GGG. The only commonality in all three systems 
is the broad feature signaling short range correlations at 
around 0.8 K in GGG and 1.0 K in GAG and GTLG. 
The lower temperature of the broad feature in GGG is 
consistent with its smaller nearest neighbor exchange in- 
teraction Ji = 0.107 K [ll|,|27| as compared to Ji = 0.126 
K for GTLG [H and Ji = 0.142 K for GAG 

The sharp features observed in GTLG and GAG are 
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states of GTLG and GAG appear to be robust against 
levels of random dilution at least as high as those in 
GGG, it seems that the type of disorder may be the cru- 
cial ingredient in the exotic behavior of GGG. In other 
words, the random excess of Gd on Ga sites, found in 
GGG, may be a much more powerful way to introduce 
random frustration and trigger a spin glass transition 
than simple dilution of the magnetic moments. This work 
shows that the response of a HFM to random disorder 
can be highly nontrivial, and a function of the nature 
of that disorder in concurrence with the nature of the 
interactions present. 
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FIG. 3: (color online) Comparison of the specific heat of GGG 
(blue squares) , GTLG (green triangles), GAG (red diamonds) 
and a sample of 2% diluted GTLG (violet x's). 



in all likelihood signatures of transitions to long range 
order. In fact the majority of insulating rare-earth gar- 
nets studied I28l430l exhibit transitions to a magnetically 
ordered state [3l|. An exponential drop in the specific 
heat below the transitions, indicative of gapped spin- 
wave excitations, also provides strong evidence of LRO. 
Such behavior is reminiscent of the Gd pyrochlore mate- 
rial Gd2Sn207 which also shows a sharp first-order phase 
transition to a ground state exhibiting static magnetic or- 
der and well-defined excitations as seen b y n eutron scat- 



tering [33 and specific heat experiments [19|. 

The materials GTLG and GAG clearly do not share 
the same glassy physics as GGG since spin glasses are 
universally found to not exhibit a sharp peak in C [ssj . 
These results prove that the glassy physics of GGG is not 
a ubiquitous property of Gd garnets and we are therefore 
left with two possible conclusions. First, a very precise 
tuning of the parameters of the Hamiltonian may be re- 
quired in order to produce the necessary competition be- 
tween local and long range interactions to give rise to 
a magnetic analog of the glass transition, even without 
quenched disorder. GTLG and GAG may be outside this 
narrow "window" of required parameter space. However, 
in many glass forming liquids (see for example [34| ) it is 
found that tuning interaction strengths, for instance by 
altering the substituents of molecules, without changing 
the overall symmetry of those molecules, does not tend 
to preclude the glass transition. 

A plausible alternative is that it is indeed a sensitivity 
of the system to small levels of disorder that results in a 
spin glass transition in GGG. Since the ordered ground 
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